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Abstract 

We discuss the non-local 'Variant Monte Carlo' algorithm which has been successfully em- 
ployed in the study of driven elastic strings in disordered media at the depinning threshold. 
Here we prove two theorems, which establish that the algorithm satisfies the crucial 'no- 
passing' rule and that, after some initial time, the string exclusively moves forward. The 
Variant Monte Carlo algorithm overcomes the shortcomings of local methods, as we show 
by analyzing the depinning threshold of a single-pin problem. 
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1 Introduction 

Over the past few years relevant progress has been made in the study of disordered 
elastic interfaces. One of the most intriguing problems concerns the response of 
the interface to an external force /. Two regimes are observed at zero temperature: 

(i) When / is smaller than a certain critical threshold f c , the interface is pinned; 

(ii) when the force / passes the threshold value (/ > f c ) the system undergoes the 
so-called depinning transition [1], which has been widely investigated during the 
last years [2,3]. The functional renormalization group has allowed to gain a much 
deeper understanding of this transition [4]. A number of experiments on the contact 
line of a liquid meniscus on a rough substrate have also been analyzed [5,6,7]. 
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In this context, we have introduced new algorithms which allow us to solve the de- 
pinning problem at zero temperature in finite samples. This paper discusses mathe- 
matical aspects of the Variant Monte Carlo (VMC) algorithm, which have not been 
published yet, although they were implied in past works [8,9,10]. In section 2, we 
prove the no-passing theorem and the forward-moving property for the VMC algo- 
rithm, then in section 3 we discuss in detail the single-pin problem. 

The VMC algorithm is able to detect the critical force and the critical configuration 
(i.e. the ultimate pinned configuration) of a one-dimensional string with short-range 
elasticity moving on a disordered two-dimensional lattice. We have also developed 
a continuous algorithm [11], which is remarkably useful in higher dimensions and 
allows to handle long-range interactions. 



2 Variant Monte Carlo Algorithm 



We consider a string h l = {/i*}i=o,...,L moving at times t — 0, 1, 2, ... on a spatial 
square lattice of side L, in a random potential V(i, h) with h = 0, . . . , M. The 
energy of the string is 

E(h*) = £ [V(i, hi) - fh\ + E d (\hli - KD] , (1) 
i=i 

where / is the external driving force and E el a short-range convex elastic energy. 
We assume toroidal boundary conditions with a winding term for / such that, at 
large /, the line keeps winding around the torus, lowering the energy at each time 
step. In section 3, we show that a non-local algorithm needs to be used allowing 
an arbitrary number of points i to move simultaneously by one site in any direction 
[8]. 

Following Ref. [8] we define a 'forward front' of length A; as a contiguous set of 
points i, % + 1, % + k — 1 which move together in forward direction. A 'backward 
front' is defined similarly. A front is 'unstable' if moving it lowers the energy (1). 
One has to check only ~ 2L 2 fronts to establish whether a string is pinned, ie has 
no unstable front. The 'depinning force' fd(h a ) of a string h a is the smallest non- 
negative / which destabilizes one of its forward fronts. The 'critical force' of a 
whole sample can then be defined as the largest of the depinning forces of all the 
strings in the sample 

/ c = max (2) 

The VMC algorithm simply moves a single front of minimal length k among the 
unstable forward and backward fronts. The VMC algorithm is not a valid Monte 
Carlo algorithm. However, each possible move within the VMC algorithm is also 
allowed with all the non-local algorithms, and its depinning threshold and ultimate 
pinned configuration is the same as the one of each non-local algorithm. We will 
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Fig. 1. Illustration of the two theorems. 

show that the above definition of f c is appropriate for the VMC algorithm. This im- 
plies that it is also correct for general non-local rules, even if they are not restricted 
to moving fronts only ([8]). 

We stress that these definitions and the following theorems can be easily extended 
to a general d-dimensional interface with a convex elastic energy. However, the 
VMC seems to be practically useful only for one-dimensional interfaces with short- 
range elastic energy as the total number of fronts remains polynomial. We prove the 
following theorems: 

Theorem 1 Let h a be a stable configuration. Then the VMC algorithm cannot 
reach a configuration h? with h] > hf for some i from a starting configuration 
hP with h? <h?Vi. 



This theorem is illustrated in figure l(a, b, c and d). We suppose the existence of a 
string li 1 forbidden by the theorem (see figure 16): 



/) E(h^) - E(h p ) < 
II) h] > for some %. 



(3) 



Due to the definition of the VMC dynamics the front connecting hP to K 1 ' (see 
figure lc) must be stable. Moreover, the stability of h a assures that the front con- 
necting h a to li 1 " (see figure \d) is also stable: 



/) E(h~<') - E^h 13 ) > 
//) E{h~<")-E{h a ) > 0. 



(4) 



Subtracting Eq. (31) from the equations (4) and using the expression (1) leads to 

M\h a m - 1|) - E d (\h a J) + E d (\hi + 1|) - E el (\hi\) > 0; (5) 



where, without any loss of generality, we set h^+i = ^m+i = ( see figure la). 
From ft£ < hZ we write: 



hi < h a m - 1< h a m 
hi. <hi + l< h° 



(6) 
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A convex function f(x) in [x, y] satisfies for all x±, x 2 G [x, y]: 

f(t Xl + (1 - t)x 2 ) < tf( Xl ) + (1 - t)f(x 2 ) 0<t<l. (7) 
Taking x x = hP m and x 2 = h^, from (6) we can find a t such that: 

h a m -l = thi + (l-t) lf m 

hi + l = th a m + (l-t)hi. (8) 

Using this relation to impose the convexity of the elastic energy in (5) we end up 
with a contradiction which demonstrates the theorem. 

Theorem 2 Let h a be pinned in forward direction. Then, the VMC algorithm can 
at most recede towards a string hP (Tif < hf which is itself pinned in forward 
direction. The analogous property holds for strings pinned in backward direction. 

An illustration of this theorem is displayed in figure l(e, /, g and h). We suppose 
that a configuration li 1 , forbidden by the theorem, is reached by the string: 

I)E{li<) - E{h,P) < 
II)h]>h« for some i. (9) 

The instability of the front connecting h a to h 13 (figure le) implies 

E(h p ) -E(h a ) < 0. (10) 

We may connect the string h a to the string li 1 by moving two fronts (see figure 
Ig) 1 . The forward front is stable because h a is pinned in forward direction. The 
backward front is also stable because it is smaller than the front connecting h a to 
hP (VMC dynamics). Thus we conclude: 

E(h?) - E(h a ) > 0. (11) 

Eq. 1 1 contradicts the sum of (9 T) and (10), invalidating the starting hypothesis and 
proving the theorem. 

The 'no-passing theorem' (theorem 1) assures that the VMC algorithm connects 
an arbitrary initial state with the critical string, whereas the 'forward-moving the- 
orem' (theorem 2) allows us to understand that Eq. (2) is indeed appropriate: one 
might have imagined that the elastic line which cannot advance at f c could move 
backwards and then be avoided during the subsequent forward evolution. Theorem 
2 excludes the existence of such loopholes. Finally we remark that both theorems 
transpose correctly to the lattice the analytical properties of the continuum equation 
of motion [12]. 

1 In analogy with theorem 1 we have also to consider the case shown in figure 1 h, where 
we suppose that the line starting from hP moves to h x instead of b? '. Using the convexity 
condition (7) the theorem remains valid. 
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Fig. 2. Single-pin problem with harmonic elasticity. The stationary solution of the VMC 
algorithm is close to the continuum solution. The stationary solution of the local Monte 
Carlo algorithm approaches the continuum solution only for / > /i oc = 1. 

3 Single-pin problem 

We now discuss the motion both of a continuous and a discrete elastic line. The line 
is pinned at a single point which corresponds, with periodic boundary conditions, 
to h — Iil — 0. The equation of motion for the continuous line is 

dh/dt = -dE A /dh + f. (12) 

A discrete harmonic elastic energy E e \ = \\h i+ \ — hi\ 2 then corresponds to the 
continuous energy E e i(x) = ^(dh/dx) 2 , to be integrated over x. The stationary 
solution under the indicated pinning condition is easily seen to be h(x) = \ fx(L — 
x), as shown in figure 2 for L = 10. 

We may also follow the dynamics of the discretized problem from a starting con- 
figuration hj =0 = Vi Of particular interest is the case / = 1 ± e, with e > 0. The 
VMC solutions are given in the figure. They have to be contrasted with the solu- 
tions of the local algorithm (only fronts of length 1 can be moved). For / = 1 — e, 
the starting configuration is stable under local dynamics, as any forward move of a 
single point costs an elastic energy 1, more than what is recovered through the driv- 
ing force. The VMC solution is recovered only for / > / loc l This specific 'critical 
force' / loc of the local algorithm is independent of L only for the harmonic elastic 
energy. 

The local algorithm is more pathological for stronger than harmonic elastic ener- 
gies, which have proven to be important in this context [9]. This is evident in the 
metric constraint model (see [8]). As an example (stronger than harmonic, weaker 
than metric constraint) we treat a quartic elastic energy E el = -^\h i+ i — hi\ 4 (corre- 
sponding to a continuous energy E el (x) = -^(dh/dx) 4 ). The stationary continuum 
solution of Eq.12: h(x) = |(^-) 1 / 3 [L 4 / 3 — (L — 2x) 4 / 3 ], is again recovered with the 
VMC algorithm. However, the specific critical force at which the local algorithm 
becomes equivalent to the VMC method grows with L as / loc ~ (|/) 2 / 3 L 2 / 3 , i.e. 
diverges with L. 
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We conclude that the local dynamics is inconsistent even for a single-pin problem. 
In general disordered samples, it similarly fails to describe the real dynamics of the 
continuum and is very sensitive to exceptional (local) configurations of the disorder 
potential, which can block the string even of an infinite systems, and eliminate the 
interplay between disorder and (collective) elasticity which is at the heart of the 
depinning problem. 
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